Temperature dependence of the slip length in polymer melts at attractive surfaces 
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Using Couette and Poiseuille flow, we extract the temperature dependence of the slip length, 
6, from molecular dynamics simulations of a coarse-grained polymer model in contact with an 
attractive, corrugated surface. 5 is dictated by the ratio of bulk viscosity and surface mobility. At 
weakly attractive surfaces, a lubrication layer forms, S is large and increases upon cooling. Close 
to the glass transition temperature, Tg, very large slip lengths are observed. At a more attractive 
surface, a"sticky surface layer" is build up, which gives rise to a small slip length. Upon cooling, 
5 decreases at high temperatures, passes through a minimum and grows upon approaching Tg. At 
strongly attractive surfaces, the Navier-slip condition fails to describe Couette and Poiseuille flow 
simultaneously. The simulation results are corroborated by a schematic, two-layer model suggesting 
that the observations do not depend on the details of the computational model. 



Rapid miniaturization of microfluidic devices has led 
to new questions in hydrodynamics [1] which often in- 
volve the role of boundaries. The microscopic structure 
at the solid-fluid interface dictates wettability and fric- 
tion and much effort has been directed towards tailoring 
these properties. The effect of microscopic surface prop- 
erties can be incorporated into a continuum description 
via a boundary condition (HBC) to the Navier-Stokes 
equation. Theoretical and experimental studies of fluid 
flow in confined systems have demonstrated that the no- 
slip boundary condition, which postulates that the tan- 
gential velocity at the surface vanishes, can be violated 
[21 and the fluid slips past the surface. In the context 
of microfluidics, shp is often preferred because it reduces 
the amount of shear stress or pressure drop required to 
maintain flow. Hence, higher through-put can be gen- 
erated in microfluidic devices. The amount of slip can 
be quantified by the slip length, 6. The slip boundary 
condition was formulated by Navier in 1823 as a balance 
between viscous stress and friction stress at the surface: 
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where A = rj/S denotes the friction coefficient and r/ the 
bulk viscosity of the liquid. Zh characterizes the posi- 
tion at which the hydrodynamic boundary condition is 
to be applied. Both parameters, S and Zh, characterize 
the HBC. The value of the slip length, 5, has been mea- 
sured experimentally for a variety of systems and the 
results have attracted abiding interest y . Typically, the 
scale of the slip length in simple liquids is set by the size 
of the constituents of the fluid 3]. There are, however, 
notable exceptions: For instance, topographically struc- 
tured surfaces have been designed to generate large slip, 
nano-bubbles or the formation of a lubrication layer at 
the solid-fluid interface due to surface segregation of a 
low-molecular weight component in a polydisperse melt 
can increase the slip length. Recently, a careful analy- 
sis of the rim profile of dewetting polymer films has also 
provided evidence of slip length of the order of microme- 



ters, which is significantly larger than microscopic length 
scales in the system \^. 

In this Letter, we investigate the impact of surface in- 
teractions and the concomitant effective viscosity [5|] at 
the surface on the HBC by molecular dynamics simula- 
tions of a standard, coarse-grained polymer model [6|. 
We show that (i) the stronger is the attraction between 
fluid particles and substrate the smaller is the slip length, 
(ii) slip lengths can grow large as the glass transition tem- 
perature of the fluid is approached, and (iii) for strongly 
attractive surfaces the Navier slip condition ^ fails to 
provide a HBC that simultaneously describes Poiseuille 
and Couette flow. Our simulation results are corrobo- 
rated by a phenomenological two-layer model. 

In our model, polymers are comprised of A'^ = 10 
coarse-grained segments that interact via a truncated and 
shifted Lennard-Jones potential with cut-off, Vc — 2\/2a. 
Length and energy scales are set by the parameters a 
and e of the potential. Neighboring segments are bonded 
together by a FENE (Finitely Extensible Non-linear Elas- 
tic) potential [q]. The model exhibits a glass transition 
temperature and the dynamics in the bulk and in con- 
finement has been recently reviewed [7|]. A Dissipative 
Particle Dynamics (DPD) thermostat is used to control 
temperature, T. The solid surface is modeled by two rigid 
layers of Lennard-Jones interaction centers arranged on 
a FCC lattice [Sj. By varying the strength, e^, of the 
Lennard-Jones potential between solid and fluid, we tune 
the adhesion and slip. 

In order to mimic experimental conditions, the simu- 
lations are performed at the coexistence pressure where 
the polymer melt coexists with a vapor of vanishingly 
low density. Rather than using a constant pressure al- 
gorithm we tune the distance, H, between the walls to 
attain the coexistence density. The system size is chosen 
large enough for the properties in the center of the film 
to be independent from the strength of the solid-fluid 
interaction. Periodic boundaries are applied in both, x 
and y directions and the dimensions of the simulation 
cell are L^ — 19.84 a and Ly — 19.93 a, respectively. 
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FIG. 1: (color online). Bulk viscosity as a function of tem- 
perature. The circles mark the results of the molecular dy- 
namics simulations and the solid line a fit using the power 
law predicted by mode-coupling theory. The inset presets the 
liquid- vapor co-existence density, pl, versus temperature, T. 



At low temperature the fluid consists of 20 000 coarse- 
grained segments and the distance between the walls 
varies between 51 a and 57 <t. At higher temperature 
ksT/t > 0.8, when the effect of the substrate propagates 
not that far into the film, we use 10 000 particles, and H 
varies between 30 a and 38 a. 

As a first result, we depict in Fig. [l]the temperature 
dependence of the bulk viscosity extracted from the ve- 
locity profile, v{z), of Poiseuille flow at the center of the 
film. The steady-state Navier-Stokes equation yields 
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where p^, rjB, and / denote the number density of the 
liquid, the bulk viscosity, and the volume force applied 
to the segments, respectively. The inset represents the 
liquid density, pL, at coexistence. Upon reducing tem- 
perature, the model exhibits a glass transition at Tg. The 
bulk viscosity above the glass transition temperature, Tg, 

iBiT)^\T-Tg\-^ 



3|. Fitting the sim- 
0.41 and 7= 1.57, 



is well describable by the power law, 
predicted by Mode Coupling theory 
ulation data in Fig.[T]yields ksTg/e = 
which agrees with previous studies [10|. 

Due to pronounced layering effects at the solid surface, 
the effective, near-surface viscosity, ijs, deviates from the 
bulk behavior. While viscosity is only properly defined in 
the bulk, there are several ways to estimate an effective 
local viscosity. One consists in computing the local shear 
stress and defining rjs as the ratio between the local shear 
stress and the velocity gradient. However, the gradient 
of the velocity profile is difficult to extract with sufficient 
accuracy and in previous studies an explicit form of the 
velocity profile has been assumed 



Here, we instead 
analyze the local mobility of the fluid. Assuming that 
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FIG. 2: (color online). Local diffusivity as a function of the 
strength of the substrate ta at ksT/e = 0.8. The circles 
and the squares are respectively the diffusivity near to the 
substrate and in the bulk. The straight line marks the Ein- 
stein estimate D — kBT/STrrja, for the fluid of unconnected 
segments. The right inset represents the effective viscosity 
ratio, tib/tis, as a function of ts, the straight line delineates 
the formation of a lubrication layer and a sticky layer at the 
solid-fluid interface. The left inset depicts density proflles, 
and the vertical lines deflne the bulk and surface regions used 
in the calculations. 



the mobility can be estimated by the Einstein relation 
D oc 1/77, we obtain the estimate r]B/r]S ex Ds/Db |la |. 
The local diffusivity is obtained from the mean 
square displacement parallel to the surfaces as Da = 
limt^oo(Ar^ )c/4t where only particles that remain in 
region a for the entire time of the calculation contri- 
bution to the average [I4I. We have computed the lo- 
cal diffusivity at temperature /csT ~ G.8e for eg = 
0.4, 0.5, 0.6, 0.8, 1.0 e. The boundaries for the bulk region 
are chosen such that the bulk properties are independent 
of the substrate strength. The fluid properties gradually 
change as a function of the distance from the solid-liquid 
interface. Our choice of the width of the boundary region 
is a compromise: It is wide enough such that segments 
remain sufficiently long in the boundary region to deter- 
mine the diffusivity, and it is narrow enough in order for 
the properties to be dominated by surface effects. As we 
decrease eg, the near-surface mobility becomes larger as 
shown in Fig. [2] and the average time a particle stays in 
the surface region decreases. In the bulk region, the dif- 
fusivity is independent from the solid-fluid attraction, es- 
In the inset, we depict the ratio Ds/Db and observe that 
substrate strengths larger than approximately 65 — 0.5 e 
result in a ratio ris/ris < 1, i.e., a "sticky surface layer" is 
formed. For eg < 0.5 e, the surface mobility is enhanced 
and a lubrication layer forms at the solid-fluid interface. 

In order to measure the slip length, 6, and position, 
Zh, of the HBC, we simultaneously compute Couette and 
Poiseuille flow profiles [8|, [l3| for the different solid-fluid 



interaction strengths as a function of temperature. The 
shear rate employed is small enough such that 6 does not 
depend on the strength of the flow IJ] . To simulate pla- 
nar shear flow (Couette) the surfaces are moved at con- 
stant velocity, vs- We typically employed vs = 0.2a/T 
where r denotes the reduced Lennard- Jones time unit. 
Poiseuille flow is generated by applying a force, / on 
all particles. In our simulations, the volume force varies 
between / = 0.001 — O.OOSe/cr. Complementary informa- 
tion is obtained by calculating the friction coefficient. A, 
from equilibrium molecular dynamics by integrating the 
transverse force auto-correlation function 15| . Previous 
studies for our model have shown that this procedure 
gives rise to consistent results at high temperatures [8[. 
The velocity in the bulk region at the center of the 
film, z = 0, are fitted by linear or parabolic profiles 
as predicted by macroscopic hydrodynamics for Couette 
and Poiseuille flow. Let zc and zp denote the positions, 
where these linear and parabolic proflles extrapolate to 
vanishing velocity, then the slip length and the position 
of the hydrodynamic boundary condition are given by: 



z^ — Zp and z/i = zc + 5 



(3) 



The results for the slip length are depicted in Fig. [3] and 
the distance to the wall of the hydrodynamic boundary 
z/i in the inset. We see that the temperature dependence 
of S is highly sensitive to temperature and solid-fluid in- 
teraction. First, we observe that independently from eg 
the slip length diverges as we approach the glass tran- 
sition temperature, Tg, because the fluid eventually be- 
haves like a solid. Already at ksT/e = 0.5, i.e., about 
20% above the glass transition temperature of our model, 
5 has increased by an order of magnitude compared to 
the approximately constant value at high temperature. 
This observation offers an explanation for the surpris- 
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FIG. 3: (color online). Slip length, 5, versus temperature 
for different strengths, eg, of solid-ffuid interaction. The in- 
set represents the position of the hydrodynamic boundary Zh 
versus temperature. 



ingly large slip length observed in the dewetting exper- 
iment of Fetzer et al. [J|, which were performed in the 
vicinity of the glass transition temperature. 

The results for the position, z;j, of the HBC do not vary 
significantly with the strength of the ffuid-solid interac- 
tion over the entire temperature regime. At high tem- 
perature, the position zt is close to the top of the solid 
surface, while it gradually shifts inwards as the tempera- 
ture is reduced towards Tg. This effect goes along with a 
growing distance over which the liquid structure is altered 
by the surface as can be observed from the pronounced 
packing effects in the density profile. 

While the position, z/j, does not significantly depend on 
eg, the behavior of the slip length qualitatively changes 
with the strength of the solid-fluid interaction. While 
the values of 5 corresponding to e^ = 0.4 e and e^ — 
0.5 e decrease monotonously with T , we observe a non- 
monotonous variation of the slip length for larger attrac- 
tion, es- For very strong attraction, eg = 1, there is even 
a region of intermediate temperatures where zc < zp 
and thus Eq. fSj has no solution. This marks the failure 
of the Navier slip condition to parameterize the near- 
surface flow solely by the material properties of the solid 
surface. 

In order to rationalize these simulation results and ex- 
plore whether they are universal or specific to our coarse- 
grained polymer model, we propose a schematic two-layer 
model depicted in Fig. ID Within this model, we ap- 
proximate the gradual variation of the fluid properties 
as a function of the distance from the solid surface by a 
boundary region of width, A, which is characterized by 
a surface viscosity, rjs, and the bulk with viscosity, tjb- 
Within each layer, the fluid is described by the Navier- 
Stokes equation. At the interface between the solid sur- 
face and the boundary layer, we impose a Navier slip 
condition ([T]) with a microscopic slip length, 5s- At the 
interface between the surface layer and bulk, we require 
the continuity of shear stress and velocity. 
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For planar shear flow one obtains for the linear velocity 




FIG. 4: (color online). Sketch of the Couette and Poiseuille 
flows in the schematic two-layer model 



profile in the central bulk region 
z + d{ 1 



V 



Vb 



rw 
VS 



Vs 



(5) 



where E quantifies the shear stress. If a volume force / is 
applied to the fluid, one obtains a parabolic velocity pro- 
file in the boundary layer and in the bulk. Using Eq. ^ 
and the boundary conditions ([T]) and ([4]), we obtain for 
the latter: 
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Finally, Eq. ([3]) yields the slip length 
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The first term describes the effect of the surface layer, 
the second term arises from the microscopic slip at the 
solid surface. For surfaces with a large surface mobility, 
Vb/vs > 1, a lubrication layer is formed and results in 
an enhanced slip length, S > ds, compared to the micro- 
scopic slip at the solid-fluid interface. On the other hand, 
if the solid-fluid interactions give rise to a boundary layer 
with large effective viscosity, vb/vs < li the presence of 
this sticky layer at the substrate reduces the slip length, 
S < Ss). Moreover, if 



r,B ^ 1 + 2(5s/A 
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the velocity far away from the surface cannot be de- 
scribed by the Navier-Stokes equation and a Navier slip 
boundary condition ([T]). 

This schematic model can rationalize the observations 
in our molecular dynamics simulation: (i) At high tem- 
perature, kinetic effects will dominate the behavior, thus 
Vs ~ Vb- In this case, S is equal to the microscopic 
slip length 5 fv Ss- (ii) Upon cooling the fluid, the 
bulk viscosity increases. If the solid- fluid interactions 
are weak, eg < 0.5, a lubrication layer is formed and 
the slip length increases, 5 oc (vb/vs) ^s- (iii) If the 
coupling between solid and fluid is strong, however, the 
ratio vb/vs decreases upon cooling and so does S. If the 
ratio becomes sufficiently small (see Eq. ^), as it does 
in the case eg = 1 for our model, the Navier slip con- 
dition fails. Upon approaching Tg from above, however, 
the slip length passes through a minimum and increases. 
The latter effect stems from the temperature dependence 
of the microscopic slip length, Ss = Jys/A, which diverges 
for T -> Tg. 

In conclusion, surface interactions can modify the near- 
surface mobility and thus can be exploited to control 
the hydrodynamic boundary condition. Very large slip 



lengths can be expected in the vicinity of the glass tran- 
sition of the fluid. Depending on the strength of the in- 
teraction between the solid surface and the fluid, slippage 
may be enhanced or reduced, and at strongly attractive 
surfaces the Navier slip condition may even fail to pro- 
vide an appropriate boundary condition to the Navier- 
Stokes equation with parameters that solely depend on 
the surface material. These findings of the simulations 
are corroborated by a schematic two-layer model which 
shows that the effects are universal and are not a con- 
sequence of the non-Newtonian nature of the polymer 
liquid employed in our study. 
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